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Recent calculations of heavy quark cross sections near threshold at next-to-next-to- 
leading order have found second-order corrections as large as first-order ones. We 
[ analyse long-distance contributions to the heavy quark potential in momentum and 

coordinate space and demonstrate that long-distance contributions in momentum 
space are suppressed as Aq CD /q 2 - We then show that the long-distance sensitivity 
of order Aqcd^* introduced by the Fourier transform to coordinate space cancels 
to all orders in perturbation theory with long-distance contributions to the heavy 
quark pole mass. This leads us to define a subtraction scheme - the 'potential sub- 
traction scheme' - in which large corrections to the heavy quark potential and the 
'potential-subtracted' quark mass are absent. We compute the two-loop relation 
of the potential-subtracted quark mass to the MS quark mass. We anticipate that 
threshold calculations expressed in terms of the scheme introduced here exhibit 
improved convergence properties. 



1. Motivation 



Heavy quark production near threshold through virtual photons or Z bosons is very 
sensitive to the quark mass and therefore may allow us to determine heavy quark masses 
precisely. Recently, the two-loop corrections to the colour Coulomb potential JTJ and 
to the matching relation between the relativistic and non-relativistic vector current 0, 
|3| were obtained. The two together provide the necessary input to compute heavy 
quark properties near threshold in next-to- next-to-leading order (NNLO). In this context 
'NNLO' means that all corrections to the Born cross section of order (a s /v) n a k s for any 
n and k = 0, 1, 2 are taken into account, where v is the small relative velocity of the 
two quarks in their centre-of-mass frame and a s is the strong coupling. (Additional 
logarithms of v are not written explicitly.) 

NNLO calculations have now been completed for top-anti-top production near thresh- 
old ||, |!|, for bottomonium threshold sum rules || [7] and for quarkonium energy levels 
||. In all three cases the NNLO correction is as large as the next-to-leading order (NLO) 
correction, which suggests that a perturbative treatment has already reached its limits. 
In case of ti production the NNLO correction shifts the location of the cross section peak 
position by about 1 GeV, which implies an uncertainty in m t of about 0.5 GeV, if the 
threshold cross section is used as a measurement of the top quark mass. This result is 
unexpected, in particular as the relevant physical scale is given by Cpa s {m t v ) m t /2 ~ 15- 
20 GeV at which perturbation theory should work. 

Following a different line of investigation, the Coulomb potential in momentum and 
in coordinate space is analysed in || [L0|. It is found (see also |pT| ) that the effective 



couplings defined by the two versions of the potential are related by a rapidly divergent 
series, the origin of which is a long-distance contribution of relative order Aqcd 7 "- This 
leads to large numerical differences in different, but consistent at NNLO, implementations 
of the Coulomb potential in cross section calculations. The authors of |S|, [TO] conclude, in 
agreement with the evidence from the NNLO ti and bb computations mentioned above, 
that there is a large and irreducible uncertainty that affects the threshold region. This 
would lead to rather gloomy prospects as to our ability to constrain the bottom and, 
eventually, the top quark mass. 

In this paper we show that despite these evidences perturbation theory does not (yet) 
fail and that the actual uncertainties can be smaller than those indicated by [[|, |5|, 0, [|, |1(| . 
Our main point is that, contrary to intuition, the notion of a quark pole mass, which 
has been implicit in the discussion above, is in fact inadequate for accurate calculations 
of heavy quark cross sections near threshold. The argument goes as follows: 

We first analyse (Sect. 2) the heavy quark potential V(q) in momentum space and 
find that long-distance contributions have a relative suppression (AQ CD /g) 2 . Hence the 
potential in momentum space is better behaved than the potential in coordinate space. 
Knowing that the long-distance contribution of relative order Aqcd?" to the potential in 
coordinate space enters only through the Fourier transform, we can eliminate it by re- 
stricting the Fourier transform to \q\ > fif for some factorisation scale /z/ which satisfies 
Aq CD < fif < mv. This defines a subtracted potential V(r, fif), from which large pertur- 
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bative corrections are eliminated (Sect. 3). The Schrodinger equation takes its conven- 
tional form only if the pole quark mass definition is assumed. The Schrodinger equation 
formulated with the subtracted potential contains a residual mass term Sm(iif). As a 
consequence the input parameter for threshold calculations in terms of the subtracted 
potential is not m po i e but mps (/•«/) = rn po \ c — 5m(fif). It is known that the pole mass 
also receives long-distance contributions of order Aq CD ]12|, [TJ| and it was noted already 
in that they are related to the Coulomb contribution to the self-energy. The crucial 
point is that the long distance sensitivity of order Aqcd in the coordinate space poten- 
tial or, equivalently, 5m(fif) cancels to all orders in perturbation theory with the leading 
long-distance sensitivity in the pole mass (Sect. 4). Hence the 'potential-subtracted' 
mass mpg(/iy) can be related to more conventional (long-distance insensitive) mass defi- 
nitions by a well-behaved perturbative expansion. The two-loop relation to the MS mass 
definition can be trivially obtained (Sect. 5). 

It follows from this argument that one can avoid large perturbative corrections as 
were found in the NNLO calculations mentioned above by formulating the threshold 
problem in terms of mps(/i/) and the subtracted potential V(r, ///) rather than the pole 
mass and the ordinary Coulomb potential. One can then determine mps(/x/) and relate 
it reliably to m^. The dependence on the factorisation scale fif cancels in this process. 

2. The potential in momentum space 

The static potential in coordinate space, V(r), is defined in terms of a Wilson loop 
Wc(r, T) of spatial extension r and temporal extension T with T — > oo [[14], [T^, [U|. 
In this limit Wc{r,T) ~ 6xp(— iTV(r)). The potential in momentum space, V(q), is 
the Fourier transform of V(r). (We use r = |r| and q = \q\.) One can compute 
the potential directly in momentum space from the on-shell quark-anti-quark scattering 
amplitude (divided by i) at momentum transfer q in the limit of static quarks m — ► oo 
and projected on the colour-singlet sector. In addition one has to change the sign of 
the ie-prescription of some of the anti-quark propagators in some diagrams (such as Di 
below), so that the integration over zero-components of loop momentum can always be 
done without encountering quark poles in the upper half plane. The quark poles amount 
to iterations of the potential, but do not give a contribution to the potential itself.Q 

The potential is infrared (IR) finite^ and ultraviolet (UV) finite after renormalisa- 
tion of the coupling. In this Section we ask how sensitive the Feynman integrals that 
contribute to the momentum space potential are to the IR regions of loop momentum in- 
tegrations. The reason is that these regions give rise to large corrections in higher orders 



in perturbation theory through IR renormalons (see |12|, [0|] for references). Note that 



we are not concerned with the long-distance behaviour of the potential at q ~ Aqcd ; but 

1 We remark that in the threshold expansion of Feynman integrals [[l7| the static potential and its 
generalisation are generated by integrating out soft quarks and gluons with energy and three-momentum 
of order mv together with potential gluons with energy of order mv 2 and three- momentum of order mv. 

2 We are aware of the analysis of fl6|| , which can be interpreted as a statement to the contrary. In 
our opinion, the interpretation of the divergences discussed in |lq| deserves further consideration. 
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Figure 1: One- loop corrections to the heavy quark potential. 



with the leading power corrections of form (Aqcd/q , ) A: , which correct the perturbative 
Coulomb potential when q is still large compared to Aq CD . 

Consider first the one-loop corrections to the potential (see Fig. |l|) in dimensional 
regularisation. In Feynman gauge D% is zero. The colour Cp-term cancels in the sum 
D\ + D 2 . D\ i are logarithmically IR divergent. This IR divergence is cancelled by the 
(logarithmic) scaleless integrals D^, whose only effect is to convert the IR singularity in 
Di j2 into a UV singularity, which can be absorbed into a renormalisation of a s . Hence, 
we are left with the IR finite parts of (the C^C^-part of) D 2 and D e , which yield the 



well-known one-loop correction to the potential |T5|, |16fl . Dq is given in terms of the 
gluon self-energy at off-shell external momentum q. For small loop momentum k, the 
integral scales as 

and gives rise to a contribution of order AQ CD /g 2 from the region k ~ Aqcd relative to 
the tree diagram which scales as 1/q 2 . The integral relevant for D 2 is 

k 2 (k + qf{v k) 27 (2) 

where v = (1,0) and v ■ q — 0. To find the leading infrared contribution, which is left 
over after the IR divergence is cancelled as described above, we expand the integrand 
in k around k = and around k + q = 0. The integrals in each term of the expansion 
depend only on the vector v. Hence, in a regularisation scheme that preserves Lorentz 
invariance all odd terms vanish because v-q = 0. The long-distance contribution is again 
of relative order AQ CD /g 2 . (Note that we are only concerned with IR contributions that 
are connected with the large-order behaviour of perturbative expansions.) 

This argument generalises to an arbitrary diagram. Because v ■ q = and because 
there is no other kinematic invariant linear in q, it follows from Lorentz invariance that 
the leading power correction to the potential in momentum space cannot be Aqcb/q, 
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but has to be quadratic: 



+ . . . + const x + (3) 

The implication is that the expansion of the potential in the coupling^ at s (q) diverges as 
E„r n a s (g) n+1 ~ J2 n (-a(3 ) n n\n b a s (q) n+1 with a = 1 and (3 = -(11 - 2n / /3)/(4vr) the 
first coefficient of the /3-function. A linear IR power correction would have led to a = 2 
and, hence, a more rapidly divergent perturbative relation. The parameter b remains 
undetermined by the above analysis, but does not influence the power behaviour. 



V(q) 



AnC F a s (q) 
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3. The potential in coordinate space 

Consider now the potential in coordinate space, defined by the Fourier transform 

nr) = j^e^V{q). (4) 

Note that in the Wilson loop definition of the potential in coordinate space there is a 
(divergent) constant related to the self-energy of the static sources. This r-independent 
term is usually discarded when one refers to the static potential, and it is also dis- 
carded, when V(r) is defined in terms of the Fourier integral above. In what follows the 
interpretation of this constant plays an important role. 

To see that the potential in coordinate space is more sensitive to long distances than 
the potential in momentum space, it is enough to take the tree approximation to V(q) 
and to calculate 

H Jf , (_^a,\ = _2£^ Aqod + 0(A3 2) . (5) 



\q\<^QCD 

It follows that the leading power correction is linear in Aqcd?" 

"31 10n 



= C F a s (e^/r) f q a ( e -^/r) 
r I Air 
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+ . . . + const x Aqcd^ + • • • 



(6) 

(7e is Euler's constant.) The implication is that the expansion of the coordinate space 
potential in a s (e~ lE r) diverges as J2 n r n a s (e~ lE r) n+1 ~ J2 n (~ a Po) nn ^ nba s(€-~ lE /r) n+1 
with a = 2, much faster than the expansion of the potential in momentum space. Note 
that in absolute terms the IR contribution is a constant of order Aqcd- 

The rapidly divergent behaviour of the coordinate space potential has been noted in 
ljj and is discussed in detail in || 0. What we add here is the observation that the 



linear power correction and, by implication, the rapid divergence originates only from the 
Fourier transform to coordinate space and is not present in the potential in momentum 



3 In this paper a s (fj.) denotes the coupling defined in the MS scheme. 
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space. Knowing this we can subtract the leading long-distance contribution and the 
leading divergent behaviour completely by restricting the Fourier integral to \q\ > [if 
with [if a factorisation scale which we make more precise later. The result will be called 
the 'subtracted potential' V(r, [if). The subtraction terms can be evaluated order by 
order in a s given V(q) to that order. The relevant calculation will be done in Sect. 5. 
To be precise, we write 

V(r,[if) = V(r)+25m([i f ), (7) 

where 

To subtract the leading long-distance contribution of order Aqcd, it is legitimate to 
replace e %q ' r by 1 in the Fourier transform and we use this for the definition of the 
subtraction term. We now define the 'potential subtracted' (PS) quark mass parameter 
as 

m PS ([if) = m po ie - 5m([i f ). (9) 

In terms of the subtracted potential and the PS mass the Green function is determined 
from the Schrodinger-type equation 



V 



2 



™ps(a*/) 



+ V(r,[i f )-E 



G c (f,E) = 5^(f), (10) 



and it is important that the non-relativistic energy is defined as E = y/s — 2mps(/ff), 
with s the centre-of-mass energy, as compared to the usual definition ^Js — 2m po \ c . The 
equation above is the conventional Schrodinger equation, but with a 'residual' mass 
term.Q The subtracted potential V(r,[Mf), which contains the residual mass, does not 
suffer from large loop corrections associated with the leading asymptotic behaviour of 
its perturbative expansion. 

Despite this fact we have not yet gained anything, because the large loop corrections 
have only been hidden in the contribution 8m([if) to mps(/i/). The crucial point is 
this: When m po i e is expressed in terms of a 'short-distance' mass parameter such as the 
MS mass my^ through a perturbative series, this perturbative series also has large loop 
corrections [12, The large perturbative corrections absorbed into 5m([tf) cancel 



exactly with large perturbative corrections to the pole mass. The argument is given in 
the following Section. Hence one can first determine the PS mass from the threshold 
cross section without encountering large corrections, because of the subtraction in the 
potential. One can then relate the PS mass to the MS mass, again without encountering 



4 We have replaced m po i c by mps(/ff) also in the kinetic energy term in the Schrodinger equation. 
The difference to using the pole mass is Sm/m'V 2 , i.e. of higher order in 1/m. Since other terms of 
order 1/m are neglected in the Schrodinger equation, the replacement of m po i e by mpg(/tj) is justified. 
If one chooses the PS mass as the mass definition and derives the Schrodinger equation from Feynman 
diagrams, the kinetic energy term is divided by the PS mass by construction. An additional finite 
renormalization of the kinetic energy term can be neglected in the present context. 
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large corrections related to the asymptotic behaviour of perturbative expansions. In this 
way, one can, in principle, determine the MS mass from threshold cross sections to better 
accuracy than the pole mass, the use of which is implied by the unsubtracted potential. 
Conversely, one can begin with m^, compute the PS mass and predict the threshold 
behaviour. The perturbative relation between the PS mass and the MS mass is given in 
Sect. 5. 

One may ask why we do not suggest to avoid using the coordinate space potential and 
to work with the momentum space potential directly. The reason is that the Schrodinger 
equation for the Coulomb Green function Gc{p,E) in momentum space contains the 
integration 

d 3 k ~ -, ~ 

' s V{k)G c {p-k,E), (11) 



(2tt) 3 

which contains exactly the same leading long-distance sensitivity as the Fourier integral 
(|), because e^' r may be replaced by 1 as far as the leading power in Aqcd?" is con- 
cerned, cf. (§). Hence the problem one encounters with the unsubtracted coordinate 
space potential enters in momentum space when one solves the Schrodinger equation. 

How large can fif be? We shall see below that the expansion of 5m(fif) is naturally 
expressed in terms of a s (/i/). Perturbativity hence requires /i/ > Aq CD . The scale rel- 
evant to the potential is 1/r ~ mv. Since the subtraction should affect the potential 
only at distances larger than the physical scale of the process described by the poten- 
tial we require also [if < mv. There is another way to arrive at this constraint. If 
v = y/l — 4m 2 ole /s is expanded in terms of mps(/i), one generates singular terms (as 
v — > 0) of order (5m(n f) / (mv 2 )) k . These terms are small if 5m(nf) is small compared 
to scale mv 2 of binding energies of a Coulomb system. Counting a s ~ v and using 
8m(fif) ~ one arrives again at \if < mv. 

4. Cancellation of long-distance contributions with the pole mass 

Expressing the pole mass in terms of the MS mass and Sm(fif) as a series in a s , we can 
write 



mps(fif) = m po ie - Sm(iif) = M 



n=0 



-A*/£ a » a ? +1 . (12) 



where M = "^ms(^ms)- Both series diverge as r n , s n ~ (— a/3 ) n n\n b with a = 2. We now 
show that this behaviour cancels in the difference in (0). Because this divergence arises 
from long-distance sensitivity of order Aqcd in the Feynman diagrams that contribute 
to the two series, it is enough to show that the corresponding lineaif] IR behaviour of the 
Feynman integrands cancels to all orders in perturbation theory in the difference. The 
remaining long-distance contributions to the difference are of order Aq CD /M and the 
corresponding divergent behaviour has only a = 1. This establishes that the PS mass 
can be reliably related to the MS mass. 



5 A loop integral that behaves as J d 4 k/k A for small k is called logarithmically IR sensitive, an integral 
that behaves as J ^k/k 3 linearly IR sensitive etc.. 
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Since we are only concerned with infrared behaviour we may work with unrenor- 
malised quantities which differ from MS renormalised quantities only by pure UV sub- 
tractions. Radiative corrections to the pole mass are given by the self-energy, evaluated 
at p 2 = m 2 ole : Am = m pole — m = £j^= mpole , where m is the bare mass and S° the 
unrenormalised self-energy.Q Consider the cancellation at the 1-loop order. As long 
as we are interested only in the leading behaviour at small loop momentum, we can 
approximate the one-loop contribution to E° by 

_, c 2 f A 7 M 0+#+rao)7M 2 r d 4 q 1 

1 F9s J (2vr) 4 ((p + q) 2 - ml) q 2 p 2 =™ 2 ~^ F9s J (2vr) 4 (v q + ie) {q 2 + ie) ' 1 J 

setting p = m v. At this order one need not distinguish between m po i e and m in the 
integrand. Taking the integration over go m the upper half plane, we obtain 

1 f d 3 q AitCpOis 



2 J (2vr) 3 -q 



2 



(14) 



which is exactly the leading-order contribution to 5m(fif), cf. (]§D, for small q. Thus the 
Feynman integrands of the integrals contributing to Mtq and /i/So cancel each other in 
the infrared region of small q. Note that the leading infrared behaviour can be obtained 
by replacing 

> —in S(v ■ q) (15) 

v-q + ie K ' K ' 

in (|T3|). In other words, the relevant IR behaviour is obtained from setting q = first 
and then from the small-g behaviour of the remaining three-dimensional integral. 

The denominator of an on-shell heavy quark propagator with momentum mv + I is 
v ■ I + I 2 /(2m). To demonstrate the IR cancellation in higher-loop order, we consider 
first the static approximation, in which the denominator is simplified to v ■ I and the 
gluon coupling 7 M to heavy quarks by v^. Hence the Feynman rules reduce to those 
implicit in the definition of the potential. The static approximation implies that we 
consider all loop momenta small compared to m and take the first term in an expansion 
in 1/m. We show that the leading IR contributions of order Aqcd to the pole mass in 
this approximation cancel exactly with those to the coordinate space potential. At the 
end of this Section, we address the question of what happens, when one includes further 
terms in the expansion of the heavy quark propagator in I and the region I ~ m, in 
which the propagator cannot be expanded. 

A general self-energy diagram in the static approximation can be written as (see 
Figure |]a) 

L rl A k N 1 

where the line momenta Zj are linear combinations of the loop momenta k m and S(li) 
contains no heavy quark propagators. Consider any one-particle irreducible (1PI) sub- 
graph contained entirely in S at fixed, non-zero, external momentum. Such subgraphs 

6 The self-energy is given by the one-particle irreducible two-point diagrams divided by (— i). 
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(<-) 

Figure 2: (a) Structure of an arbitrary self-energy diagram. (b,c) Some 2-loop examples. The 
double line denotes the static approximation. 

are IR finite and at most quadratically IR sensitive. It follows that any subgraph that 
can give rise to linearly IR sensitive integrals must contain at least one heavy quark 
propagator with v ■ li — > 0. Let us call a self-energy diagram (ir-)reducible if it contains 
(does not contain) a self-energy subgraph. 

Consider irreducible diagrams such as the diagram depicted in Figure |||b first. Irre- 
ducible diagrams are IR finite on mass-shell and none of the heavy quark line momenta 
coincide. The leading IR behaviour is obtained by letting one of the v ■ ij = i° go to 
zero and by neglecting Z° in the other propagators. This can be summarised as the 
substitution 

N i N N i 

n^--E(-")^-w.n ("I 

1=1 1 J=L 1=1,1=?] 1 

In a term with S(v ■ lj) change one loop integration variable to q — lj(k m ). The delta- 
function kills the go integral and sets go — in the integrand. One can interpret (|T7j) 
as cutting the diagram at any heavy quark propagator. This gives rise to N four-point 
diagrams, to be integrated with d 3 q/(2ii) 3 , and it is easy to see that the result matches 
exactly with contributions to 5m(^f). For example, the diagram of Figure |||b cancels 
with the contribution to 5m(fif) from D 2 + 2 ■ D4 in Figure [l|. The signs and z's work 
out correctly and the factor 1/2 in ([5]) comes from the fact that we have 7r in ( |TTD but 
one factor l/(27r) from the four- dimensional integration measure. 

One may be concerned about the fact that after setting g = 0, the integrals be- 
come IR divergent and that hence the leading-order IR approximation ([T7|) may not be 
sufficient for linearly IR sensitive contributions. The IR divergences are just the IR diver- 
gences in individual contributions to the potential mentioned in Sect. 2, which cancel in 
combinations such as D 2 + 2 • D 4 . Moreover, the next term in the small- loop momentum 
expansion of the integrand is quadratically IR sensitive as shown in Sect. 2, and this is 
enough to guarantee that ([T7|) is legitimate. 

One may also be concerned about the fact that application of ([H]) does not lead to 
D 2 literally, but to an integrand which differs from (0) in that the two factors ofv-k have 



8 



different ie-prescriptions. However, in Sect. 2 we have shown that linearly IR sensitive 
contributions originate only from the IR behaviour of the Fourier integral. Hence we 
should consider q small at fixed k or both q and k small simultaneously (but not small 
k at fixed q), in which case the potential pinch singularity is not a problem. 

The situation is more complicated for reducible diagrams such as the one in Figure 0c. 
For reducible diagrams some of the heavy quark line momenta coincide and cutting such 
a quark line in the sense discussed above leads to one-particle reducible contributions 
to the potential, i. e. to lower-order contributions to the potential multiplied by on- 
shell renormalisation of the external legs. (D§ in Figure [l] is an example.) Moreover, 
reducible diagrams are IR divergent when evaluated on-shell, while AM is IR finite. 
This is related to the fact that AM is not given by Ejy =mo but by S{^ = ™ ole - To make 
the IR finiteness explicit, the contributions from reducible diagrams to AM should be 
combined with contributions at the same order in perturbation theory that arise from 
expanding Ej^ 



"° in AM: 

=m p ole 



yO 

W=m pole 



BY 



rf=m 



'as 



mo 



+ 



9, V 



x 2 f) 2 Y° 

o Y ° n I 

w=m y= 



d$ 2 



■n H) 



+ 



(18) 



This expansion reproduces precisely the combinatorial structure of self-energy subgraphs 
in reducible diagrams. Combining the various terms on the level of integrands, the re- 
sulting integral is IR finite. Moreover, after this cancellation all heavy quark propagators 



have different momenta and one can again use ([171). The 'cut' diagrams then cancel again 
with diagrams to the potential. 

Let us illustrate this for the diagram of Figure 0c. According to (|T8| ) we combine the 
diagram with the product of 1-loop contributions to the second termF] on the right hand 
side of (UP . This gives the following contribution to AM: 



d A k\ d 4 k2 
(2tt) 4 (2tt) 4 

d A ki d 4 k2 



1 



k\ k\ (v ■ h) 2 v(h + k 2 ) 
1 



(2tt) 4 (2tt) 4 k\k\v-k x v (h + k 2 )v k 2 



d 3 q 
(27p 



(27T) 



4 • - • D 5 

2 



d A k 
D 1 



k\ k\ (v ■ ki) 2 v ■ k 2 



1 



k 2 q 2 (v ■ k) 2 k 2 (k + q) 2 (v ■ k)' 



(19) 



To arrive at the last two lines we have used (|17D . Since 1/2 ■ D 5 is the on-shell wave 
function renormalisation for a single external quark leg times the leading-order potential, 

7 Further terms contribute only at 3-loop order and beyond. 
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we obtain the desired cancellation with IR contributions to the potential. The example 
and the structure of make it transparent how the IR cancellation for reducible 
diagrams extends to all orders. 

Let us return to the validity of the static approximation. Consider first the contri- 
butions to the pole mass from the region of loop momentum where all momenta are 
small compared to m. All heavy quark propagators can be expanded about the static 
limit. The corrections to the static approximation are suppressed by at least one power 
of loop momentum divided by m. This implies a suppression of long-distance sensitivity 
by a factor of Aqcd/^ relative to the leading term. Since the leading term is already 
linear in Aqcd, we conclude that if the heavy quark line momentum is small, it is suf- 
ficient to keep only the leading term in the expansion of the propagator. Consider now 
the contributions from the region of loop momentum where some loop momenta are of 
order m and others (at least one) are small compared to m. The hard subgraphs with 
loop momentum of order m reduce to local interactions of form (li/m) k with respect to 
the small loop momenta U ~ Aqcd- One obtains contributions suppressed by powers 
of Aqcd/^ unless k < 1. Hence we need to consider only hard self-energy and vertex 
subgraphs. The effect of these subgraphs is to renormalise the coefficients of the ip^ipA 
and ip'doi/j interaction terms in the non-relativistic effective Lagrangian, from which the 
potential is derived. In the standard normalisation of the non-relativistic Lagrangian the 
coefficients of these operators are 1 to all orders in perturbation theory. It follows that 
the hard subgraphs have no effect on the IR cancellation, or, in other words, they are 
implicitly taken into account through the coefficient functions of the interaction terms 
in the non-relativistic Lagrangian which enter the calculation of the Coulomb potential. 

5. Relation to the MS mass definition 

It is straightforward to work out the mass subtraction 8m{fif) from known results on 
the potential in momentum space: 



v{q) = ¥ — 



47T \ 4:71 



(20) 



with ai as in (|) and a 2 = 634.402 - 66.3542^/ + 1.246nJ 0. From the definition (|) 
one obtains 



C F ot s {ii) 




i+ ^(--43- 2 i) (21) 



^/ _ o\ i ( 1„2 /4 _ A ,„ /4 



a 2 - {2 ai b + h} In ^ - 2 + b 2 In 2 ^ - 4 In ^ + 



o I 1 U I 9 9 



where b = -4tt/3 = 11 - 2n//3 and b x = -(4tt) 2 /3i = 102 - 38n//3. Note that the 
logarithms disappear when the coupling is normalised at the scale fif, which follows from 
the fact that the potential is physical and independent of fi. At /if — 1-1.5 GeV, a typical 
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scale relevant for bottom quarks, the third-order term already exceeds the second-order 
term. This is not a point of concern as the series expansion of 8m{fjLf) is expected 
to behave badly and we are interested only in the combinations V(r) + 2<5m(/i/) and 
m poie — 8m(fif), both of which have better behaved series expansions. The subtracted 
Coulomb potential V(r) + 28m(nf) at leading order in a s is given by 

C F a s (/i r ) f 2 a,(|i/) 



V(r) 



TT a s (n r 



■fifT 



(22) 



with fi r = e~ lE jr. To see the numerical effect of the subtraction, we choose the values 
r = 1/(20 GeV) and /!/ = 3 GeV (rif — 5) as would be appropriate to tt production and 
compare the subtracted and unsubtracted Coulomb potential. The result, including the 
known higher-order corrections, is 

2 



V(r,n,) 
as compared to 
V(r) 



C F a s (jj, r 



0.86 + 0.16 



71 



+ 13.64 



+ 



71 



■ CFa ' M (l + 1.19^ + 32.93 ( 
r \ n \ 



+ . . 



7T 



(23) 



(24) 



The convergence of the series is improved and the strength of the potential is reduced. For 
bottom systems one observes a similar effect, although the requirement that /i/ > Aqcd 
does not allow us to choose /i/r as small as we would like. 

Since the relation of the pole mass to the MS mass is known only to second order 
I8| , we can only make use of 8m{^f) to second order to express the potential-subtracted 
mass mps(/i/) = m po \ c — 5m(fif) in terms of M = ^ms(^ms)- ^ ne resm t i s 

4a a (M) r _ 
3vr [ M_ 



mps(tif) =11 + 



+ 



'a s {M) 



2 r 



7T 



K 



3M 



ai - b 



+ ... 



(25) 



H f /M 


' 1 


^PS 
'2 





1.33 


13.44 


- 1.04n/ 


1/25 


1.28 


12.07 


- 0.95n/ 


1/5 


1.07 


8.927 


- 0.73n/ 


1/4 


1 


8.207 


- 0.68ri/ 


1/3 


0.89 


7.165 


- 0.61n/ 



Table 1: First and second-order (in a s {M) /tt) coefficients in the relation between the PS 



mass and the MS mass (25). For \ij = the PS mass coincides with the pole mass. 
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where K = 13.44 — 1.04n/ fromP] [IT3 and a± = 10.33 — 1.1 In/. Numerical values for the 
first and second-order coefficients for various values of Hf/M are given in Table[l|. For 
small values of /i//M as relevant to ti production the series is not very different from the 
series for the pole mass, reflecting the fact that at this order both series still converge 
well. In absolute terms the difference between the pole and PS mass may still amount 
to several hundred MeV, which is significant close to threshold. In fact the effect of 
the subtraction is far from small on the potential even for ti production as illustrated 
by (p3|) above, because /i/ has to be compared to the scale m t v/2 ~ m t a s /2 in case of 
the potential. For ratios of /i//M that may be contemplated for bottom quarks near 
threshold the second-order coefficient in ( pip is already considerably smaller than the 
one in the relation of m po i e to M and the convergence of the series expansion is im- 
proved. This lends support to the hypothesis that the IR cancellations between m po i e 
and Sm(nf) occur not only asymptotically but already at second-order. Eq. ( p3|) and 
(p5|) taken together suggest that threshold calculations at NNLO formulated in terms of 
the subtracted potential and the PS mass exhibit reduced NNLO corrections and that 
the PS mass can indeed be reliably related to the MS mass. 



6. Conclusion 

In this paper we proposed that perturbative calculations of heavy quark properties near 
threshold should not use the pole mass but a subtracted mass together with a sub- 
tracted potential. This eliminates one source of large corrections in perturbation theory, 
related to small momentum contributions, although we cannot exclude large corrections 
due to other reasons. The numerical estimates presented above suggest, however, that 
the convergence is indeed improved. We emphasise that the potential-subtracted mass 
is gauge- invariant, because it is derived from the pole mass and an integral over the 
momentum space potential, both of which are gauge- invariant. 

The crucial point is that heavy quark cross sections near threshold are in fact less 
sensitive to long distances than the quark pole mass parameter. This follows from the 
observation that the Coulomb potential in momentum space is less sensitive to long dis- 
tances than the potential in coordinate space and that the large perturbative corrections 
to the potential in coordinate space cancel to all orders in perturbation theory with those 
to the pole mass, because of an exact cancellation of the small momentum behaviour 
of the respective Feynman integrals. That the pole mass is not relevant for physical 
quantities involving top quarks is quite obvious, because the width of order 1.5 GeV 



provides an intrinsic cut-off for long-distance effects [I9|. In particular, the location of 
the resonance-like peak in the production cross section is not a direct measure of the top 
quark pole mass despite the fact that top quarks do not hadronise. It is however less 
obvious that the pole mass is not even relevant for (quasi-) stable quarks near threshold. 

Making use of the results of [0 we derived the mass subtraction term to order 
and the relation between the potential-subtracted (PS) mass and the MS mass to order 

8 Note that mj^g(m po i e ) is used in |l8| and the different normalisation point for the quark mass affects 



the second order coefficient in the relation to 



^polc ■ 
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a 2 s . These relations provide the link with other physical quantities involving heavy 
quarks. In particular, they allow us to determine directly the bottom and top MS 
masses from bottomonium sum rules and the ti cross section, thus obviating large NNLO 
corrections that appear when these quantities are expressed in terms of pole masses 
B H> S 0- One may hope that this leads to more accurate quark mass determinations 
than for the pole masses, whose accuracy is limited to order Aqcd by long-distance effects 
fi~2| , |T3|1 independent of the process utilised to determine them. We will report on these 
applications in a forthcoming publication. 

One cannot use the MS masses themselves for threshold problems, because they differ 
from the pole masses by an amount of order ma s . This causes singular terms of order 
(a s /v 2 ) k to appear in perturbative expansions. The all-order resummation of these terms 
leads one back to the pole mass. It is necessary to introduce a factorisation scale [if and 
to choose a mass definition (the PS mass) that differs from the pole mass by an amount 
smaller than the typical energies of a Coulomb system, while at the same time not being 
too sensitive to confinement effects. This leads to a linear dependence on the subtraction 
scale (if. The use of a heavy quark mass with a linear factorisation scale dependence 
has been repeatedly advocated by Bigi et al. (see [[L3[] and the review ||20|| ). In [21] a 



mass subtraction term A(/zy) (the analogue of our 5m(/j,f)) is derived to order a s from 
certain integrals over the spectral densities of heavy-light quark currents. This mass 
subtraction differs from (§) already at order a s . This does not imply an inconsistency, 
since the necessary requirement is only that the long-distance sensitive regions cancel 
asymptotically in large orders. On the other hand, it seems to us that a subtraction 
based on the heavy quark Coulomb potential is most natural (and technically simplest) 
not only for threshold problems involving two heavy quarks, since, as observed in [ 13|| , 



the leading long-distance sensitive contribution to the pole mass is in fact conceptually 
related to the Coulomb interaction. 



Note added: After this paper was completed, we received Ref. |]22| , which addresses re- 
lated questions. The authors also note that linear IR sensitivity cancels in 2m po i e + V(r) 
and demonstrate this at the 1-loop order. 

Acknowledgements. I thank G. Buchalla, A. Signer and V.A. Smirnov for comments on 
the manuscript and N. Uraltsev for correspondence. 
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